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The scalar-isoscalar field of an effective chiral Lagrangian transforms differently in either 
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linear or non-linear frameworks: in the former case it is the counterpart of the pion whereas 
in the latter it is chiral invariant on its own. We compare the predictions from these two 
models for nucleon interactions and find results which are identical for two-body and rather 
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CN ' different for three-body potentials. Some qualitative features of three-body interactions are 
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discussed. 



O-*. I. INTRODUCTION 



Since the fundamental work of Yukawa, the meson exchange theory has been successfully 
applied to the description of the nucleon-nucleon force. Nowadays it is well established that 
the one pion exchange potential (OPEP) is dominant at large distances, yielding a strong 
force with spin and isospin dependences. 

The next layer of the interaction is associated with the exchange of two pions. This 
system has the lightest mass beyond the OPEP and involves an intermediate pion-nucleon 
(irN) amplitude in a kinematical region which is not directly accessible to experiment. 
Therefore a proper theoretical treatment of this component of the potential requires the use 
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of information based on dispersion relations, as pointed out long ago by Cottinghan and 
VinMau The implementation of this idea led to the Paris potential 0, which is rather 
successful in describing experimental data. 

An important feature of the nN amplitude is chiral symmetry, as emphasized more 
than 25 years ago by Brown and Durso J3|. Recently the application of chiral symmetry 
to the NN interaction has deserved much attention, initially in the restricted framework of 
pion and nucleon degrees of freedom 0], and nowadays it is fair to say that the theoretical 
formulation of this sector of the interaction is free of ambiguities. However, a Lagrangian 
containing only pions and nucleons cannot reproduce low-energy 7rN data and hence this 
very reliable part of the model has to be complemented mj. The combination of chiral 
symmetry and experimental information about the intermediate irN amplitude, based on 
dispersion relations, was brought into this problem in the last two years |J , with a successful 
description of asymptotic NN scattering data [0,|J. 

In contrast with the OPEP, the two-pion exchange potential (TPEP) produces an at- 
tractive interaction that depends little on spin and isospin, because it is dominated by the 
exchange of a scalar-isoscalar system. In many phenomenological potentials, the actual two- 
pion dynamics is simulated by the exchange of an effective scalar-isoscalar meson, with a 
mass around 550 MeV. In the framework of chiral symmetry, this meson is usually identified 
with the counterpart of the pion in the linear sigma-model and one obtains a prediction for 
its coupling constant to the nucleon. The existence of a scalar-isoscalar meson is contro- 
versial, the main candidate being the f (400-1200), that may be present in tctt scattering 
||. In the case of NN scattering, the TPEP is quite well accounted for by non-linear chiral 
dynamics, constrained by experimental information, and the scalar-isoscalar meson is un- 
necessary. According to the unwritten law of quantum mechanics, stating that processes 
which are not forbidden are compulsory, the exchange of two pions must be considered in 
any realistic description of the NN interaction and the use of the sigma field to simulate the 
actual TPEP gives rise to shortcomings. In particular, a scalar field with mass m s yields a 
central NN potential proportional to e~ msr /r, whereas the spatial dependence of TPEP is 
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closer to (e _Mr /r) , /j, being the pion mass Moreover, the TPEP is proportional to g 4 , 
where g is the ttN coupling constant, whereas the exchange of a sigma is proportional to 
g 2 . In spite of all these problems, the use of an effective scalar field may be useful in prob- 
lems where simplicity is more important than precision. As far as the former is concerned 
an effective scalar field allows calculations at tree level, whereas the exchange of two pions 
involves loop integrations. 

If one is willing to use an effective scalar field in a calculation, there are two possibilities 
at hand. The first one consists in employing the usual sigma field of the linear model, 
which is the chiral partner of the pion. The other possibility is to use a scalar field in the 
framework of non-linear Lagrangians, which is chiral invariant and appears naturally when 
the non- linear fields are obtained from the linear ones ||10|| . These two scalars fields couple 
differently to pions and nucleons and hence lead to predictions which do not overlap for 
some specific processes. The main purpose of this work is to explore these predictions in 
the case of three-body forces. Our presentation is divided as follows: in sect. 2 we compare 
linear and non-linear Lagrangians and in sect 3 we motivate the physical predictions of the 
latter. In sect. 4 we derive the kernel of the three-body force, which is fully developed in 
sect. 5. Finally, in sect 6 we discuss qualitatively some features of the force associated with 
sigmas and chiral scalars. 



II. LAGRANGIANS 

In this section we introduce the Lagrangians describing both the sigma and the chiral 
scalar. In the framework of linear dynamics the field a is the chiral partner of the pion field 
7r since, for an axial transformation 5 A , we have 5 a tz — > a and 5 A a — > —it. In the case of 
non-linear dynamics, the pion field is represented by <f> an d we also consider a new field S, 
which corresponds to a generalization of the field a', introduced long ago by Weinberg [JTOj . 



The corresponding axial transformation are 5 A (f> =F((p 2 ), where F(<p 2 ) is a function of <fi 2 
and 5 A S = 0. The last transformation implies that S is a chiral scalar. 
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The linear Lagrangian C a has the usual form 



+Ni fiN - gN (a + it ■ tt7 5 ) N 



7T 



\Un 2 a + U ( 



O" + 7T 



(1) 



where N is the nucleon field that transforms linearly, f n is de pion decay constant, fi is 
the pion mass, g is the ttN coupling constant and the term U (a 2 + 7r 2 ) represents the self 
interactions of the mesonic fields. In the framework of the linear sigma model the scalar 
fluctuations are associated with a field e, related to a by 



cr = U + e. 



Replacing this into the Lagrangian ([J), we obtain 

C« = \ {d^e - m 2 e e 2 ) + ~ (d,7T ■ d»n - /i 2 7r 2 ) + U ((/ w + e) 2 + tt 2 
+Ni jdN-gNtU + e + iT-n 75) N. 



(2) 
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In the non-linear approach, the Lagrangian £5 is written as 



(d^Sd^S - m 2 s S 2 ) - U (S) 



+ 



1 



2 y 



+ £n, (4) 



where / corresponds to the function / = J f% — (f> 2 . Formally, U(S) describes the self 
interactions of the scalar field and Cn represents both the nucleon sector and its interactions 
with bosonic fields. This part of Lagrangian may be cast in many different forms, two of 
which are widely employed in the literature. In one of them, the pion-nucleon coupling is 
pseudo- vector (PV) and, in the other, it is pseudo scalar (PS). For PV coupling, one has 



c pv 



fj 



N - ^^i) - mi/ii/i + —ipj^Tip ■ D^cf) - g s Sipip 
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(5) 



where ip is the nucleon field that transforms non-linearly, m is its mass and g s represents 
the coupling of the nucleon to the scalar. In this expression the pion and nucleon covariant 
derivatives are given by |TT| 
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This Lagrangian has been used recently in the study of tcN form factors in constituent 



quarks models |T2fl and it is worth noting that its last term describes a chiral invariant 
coupling between the scalar and the nucleon. 
In the case of PS coupling, one has 



c ps 



9s 



v - NijdN - 9 N(f + ir- 75 ) N - fSN (f + it ■ 75 ) N, 
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where iV is the nucleon field that transforms linearly. The last term of this expression has 
the same meaning as the corresponding one in eq. (^). 

On general grounds one knows that, in the framework of chiral symmetry, results should 



not depend on the choice of Cn ]T3|,|T4|. The equivalence between the PS and PV La- 



grangians was verified explicitly in the case of TPEP [El and of pion-nucleon interactions 



in constituent quark models ||12|| . Because of this equivalence and of the similarity with 
the linear Lagrangian, our discussions are set in PS case. Thus our complete non-linear 
Lagrangian is written as 

£s S = ~ (d^S - m 2 s S 2 ) + i (d„<(> ■ + dpfPf) + Ufi 2 f + U (S) 

+NijdN -(g + ^S^j N(f + ir- <f> l5 ) N. (9) 

The interaction terms in eqs @ and (0) have the same scalar-nucleon and pion-nucleon 
vertices depicted in fig. la. However, the non-linear approach also yields extra vertices, 
representing seagull scalar-pion-nucleon and pion-pion-nucleon interactions, among others, 
displayed in fig. lb. These results indicate that, at tree level, both Lagrangians produce 
the same NN potential. On the other hand the extra scalar-pion-nucleon vertex represents 
a genuine difference between the two approaches and has consequences in processes such as 
pion absorption by a two-nucleon system or three body forces. 



III. EFFECTIVE SEAGULL 



The purpose of this section is to discuss the meaning of the effective seagull interaction 
in terms of more basic processes. In order to do this, we note that, in the framework of 
non-linear Lagrangians, the TPEP in the pure nucleon sector is given by the five diagrams 
of fig. 2a ||, which we may want to associate with an effective scalar exchange. When an 
external pion is attached to these processes, we have the three possibilities indicated in fig. 
2b, the last one corresponding to an effective seagull, whose dynamical meaning is given in 
fig 2c. 

The nine diagrams associated with the effective seagull interaction can be understood as 
arising from the product of the amplitudes T ba and P fcca ; where the former describes the pion- 
nucleon scattering Tr a (k) N (p) — > n b {k') N {p') and latter represents the contribution from 
pion production, 7r a {k) N {p) — > n b {k') tt c (q) N (p'), as indicated in fig 3. The composite 
amplitude is denoted by A and given by 
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where Q = \ (k — k') . 

Using the ttN vertices obtained from eq. (P), we have 
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P ca = ig {(S ab T c - S bc r a - 5 ac T b + ie bac ) 
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The part of A corresponding to a scalar-isoscalar exchange is associated with the factor 
S ab m the amplitude and hence it is proportional to 

^ _ h 



5 ab P bac 



ig 3 r c u 



(p' + k') 2 — m 2 ^ 5 (p + k) 2 — m 2 
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In order to identify the main contribution to this amplitude, we go to the soft pion limit, 
by using k — k' — (k , 0) and then take the limit k — > 0. We thus obtain 

Sab p bca soft _ (13) 

which has the same structure as the effective scalar-nucleon seagull vertex predicted by 
the non- linear Lagrangian given by (H). 



IV. EFFECTIVE irNN VERTEX 

As far as NN interactions are concerned both the linear and non-linear Lagrangians, 
given by eqs. @ and (|9]), yield potentials at tree level, which are due to the exchanges of 
one pion and one scalar meson associated with the diagrams of fig. la. In principle the 
scalar coupling constants predicted by these Lagrangians could be different. However, we 
are interested only in dynamical distinctions between the sigma and the chiral scalar and, 
from now on, we set g s = g in order to force both models to simulate exactly the same 
two-body potential. 

On the other hand, these models produce different predictions for the effective 7rNN 
vertex and in this section we consider the form of this interaction. The basic ingredients are 
the amplitudes T a and T s , that describe the processes ttN — > aN and ttN — > SN respectively. 
The amplitude T a is determined by just the first two diagrams of fig. 4a whereas T s also 
includes the seagull term. As the calculations of T a and T s are quite similar, we denote 
both amplitudes by T and identify the seagull contribution by a parameter A, such that 
T a = T (A = 0) and T s = T (A = 1). At tree level, T is given by the diagrams depicted in fig 
4a, and we have 



T (A) = -%g 2 r c u fry 



b d + m £> x + m t A 



75 + 75^ o + —rl5 



pj -m 2 p 2 x - m l gf n 

with pd = p + k, p x — p' — k. Using Dirac equation, we rewrite ( |i~4"D as 



u(p). (14) 



T (A) = -tg 2 r c u (p 



k k A 
A h 



p d — m 2 pi — m 2 gfn 



7 5 m (p) . (15) 
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It is worth pointing out that, as shown in | 12| , this result does not depend on our choice 
of the PS coupling in non-linear Lagrangian and it would be the same if the PV scheme 
were adopted. 

This amplitude contains positive frequency states, which do not enter in the construction 
of the proper irNN kernel. In order to isolate these contributions, we write the nucleon 
propagator as 



fa + m 1 
p 2 — in 2 IE 



p°-E 



(p)S'(p) 



p° + E 



^v s (-p)v s (-p) 



(16) 



where E = \fp 2 + m 2 . Thus 



p + m 


1 


p2 _ m 2 





fa + m _ 
p°-E 7 



(17) 



and the positive energy contribution is 

T (+) = -ig 2 T c u (p' 
l 



2E d (p° d -E d ) 2E x (pl-E x ) 



2Eh. 2E T , 
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75« (p) 



(18) 



with pi = (Ei,pi) for i = d,x. 

Subtracting ( |18D from (|l5|) we get T, the irreducible positive frequency amplitude 



f (A) =-ig 2 r c u{p') 
( 1 



+ 



1 



2E d (p° + E d ) 2E X (jpl + E x ) 

n A 



\2E d 2E X 
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15U (p) 
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This allows the proper kernel for the process nNjNk as shown in fig 4b, to be written as 



K{\) =-tg 3 ^T c u(p>) 
( 1 1 



2E d (p° d + E d ) 2E x (p° x + E x ) 
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\2E d 2E, 
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o A 
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q 2 — ml 



u (p') u (p)] 
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where q = p' 2 — P2- 



S 



V. THREE BODY FORCE 



In this section we derive the three-nucleon force due to the simultaneous exchanges of one 
pion and one effective scalar meson and the basic diagram is shown in figure 5. It produces 
the following three-nucleon amplitude 

1 



k 2 -/i 2 



x |« (p') 
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In order to obtain the non-relativistic limit of T ljk (\), denoted by t ljk (\), we use 
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u (p') u (p) -> /, 
u (p ) 75M (p) -> A CT ■ (p - p') , 

u (p') 7°75^ (p) -> A<r • (p + p') , 

2m 

« (P') /^75^ (p) -»• -<r • k, 
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P-/1 2 " fc 2 + /i 2 ' 
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g 2 — m 2 



q 2 + m 2 ' 



where the arrows indicate that the normalization of the spinors were also changed. Using 
these results and keeping only terms of the order p/m, we get the non-relativistic amplitude 



t ijk (X) 



9 4 (0 . T U) . k 1 



(2m)- 

X<T (J) • 



k 2 + /i 2 



1 



q 2 + m 2 



(22) 



where k = p; — p'^and q = p'^, — pt- The full non-relativistic amplitude t^N is given by the 
permutation over all indices ijk. 

In momentum space the three body potential W is defined by 



(P1P2P3 \ W \ P1P2P3} = -O") 5 (p'i +P2+P3-P1-P2- Pa) tw- 



(23) 

We apply Fourier transform to ( p3|) in order to obtain the potential in configuration 
space and we have 



( r i r 2 r 3 \ w \ r i r 2r 3 } = S (r[ - ri) 5 (r' 2 - r 2 ) S (r' 3 - r 3 ) W. 

and the component W^ k of the local potential W is given by 

W ijk = 9 4 ^ m s T (i) . T (j) 
(47r) 2 (2m) 2 m 

l ~jf) U K r *i) ■ V ji<r® ■ V jt U (/ir 



(24) 
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where U(/ir) 



exp(— /xr) 



+ • V (firji) a {j) ■ V kj U (m s r kj ] 

and r kj = r k - r,-. 



(25) 



We assume the approximate identity jjh — 1 and, similarly to the procedure used in the 



two pion exchange three- nucleon force fllS] , we regularize the mesonic exchanges with the 
following dipole form factor 



'A 2 -/i 2 \ 2 



v a 2 -p; ' 

where A is a cut off parameter. 

The regularization changes the Yukawa-type function into 



(26) 
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jjr \i Ar 2 A I /i 2 
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and its derivatives are given by 

dU (/ir, A) 
<9r Q 



H—Ui {fir, A) , 



<9 2 t7 (//r, A) /1 2 

— . T [^([/( F) A)-G(r)) 
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where 
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Substituting these results into the potential, we have 



W i jk = 1 (gtV _HL_ r ^i T d) . T u) 

3 \2mJ (An) 2 \x 

x {(1 - A) [<r« • <t ( % (m s r fej -, A) ([/ (fir Jt , A) - G (r jh A)) 
+ Sij (fji, Vji) U (m s r kj ,A) U 2 {iJ>r jh A)] 



+A 



(0 . <x (J) 



xC/i (nrji,A) Ui (m s r kj , A)} , 



(33) 



where 



(fji, ffcj) = 3cr W • f jicr (j) • r fej - • f fei <x W • <r W . 



(34) 



VI. DISCUSSION 

Initially, we discuss the role of chiral symmetry in the results of the previous section. 
Inspecting eq.(25), one notes that all the terms of the potential contain two gradients, 
reflecting the fact that they come from eq.(22), which is a uniform second order polynomial 
in meson momenta, as expected from a calculation based on chiral symmetry. This feature 
of the problem is independent of A. On the other hand, the detailed form of the potential is 
quite sensitive to this parameter. The value A = corresponds to a scalar field, denoted by 

11 



a, which is the chiral partner of the pion. The value A = 1, in turn, indicates processes based 
on a chiral invariant field S. It is worth pointing out that results with A = 1 they 
should, independent of the representation adopted for the pion field and, in particular, of 
the use of either PS or PV pion-nucleon couplings. So the parameter A describes dynamical 
processes which are genuinely different. 

At tree level, these effects cannot be distinguished in NN interactions, but they manifest 
themselves in the reactions NN — > 7rNN, 7rd^->NN, in axial form factors of nuclei and in 
three-body forces. In the present work we have considered only the last kind of application. 

Our results with A = 1 coincide with that presented by Coon, Pena and Riska fl6| . These 



authors also obtained a contact interaction, employing a PV pion-nucleon coupling and a 
scalar-meson coupled to nucleons. In the framework of chiral symmetry, this combination 
means that they have tacitly used a Lagrangian equivalent of our eq.(5) and hence their a 
field corresponds to our S. 

In this work we are concerned mainly with the differences between the cases A = and 
A = 1 for three-body forces. The full exploration of this aspect of the problem would require 
a precise calculation of trinucleon observables, especially the binding energy. In order to 
produce a preliminary qualitative indication of the role of the force considered in this work, 
we evaluate the expectation value (S\ W \ S), where \S) represents the principal S-state of 
the trinucleon, which is the basic component of its ground state wave function. 



This state is written as 15 



\S) = S(r ji )Tl™(a) (35) 
where S (r,j) is the spatial component and Tl™ (a) is the antisymmetric spin-isospin wave 

2* 

function with z-components m and t respectively. 

The action of the tensor operators Sij over \S) results in states with orbital angular 
momentum different from zero. Using 



1 it 



f T (i) . T (i) a {i) . a (j) T > ( a ) = _ 3 (36) 

2* 
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we obtain 



r!r (a] 

2 1 



' W^Tl™ (a) = -4 (it 

2* 



\2mJ (4tt) 2 A* 
x {(1 - A) U (m s r kj , A) [U (//t> A) - G (/xr ii; A)] 



+ cos6j\—Ui (firji,A)Ui (m s r kj ,A) 



(37) 



where cos 0j = fji ■ r k j. Using the results of [15[ and neglecting the very short range function 
G (nrji, A) [17|], we have 



rf™ (a) ^ fc r?7 (a) = - {(1 - A) CM (m s r kj , A) [/ (/ir ii; A) 

2 l J 2 l 

+ AG S cos ^C/i (/i^i, A) U x (m s r kj , A)} 



(38) 



where C a and G5 are the strength coefficients of sigma and chiral scalar respectively, given 
by 



C a — 4 



G s = 4 



gay 



1 m s 



2m/ (4tt) /i 
4 1 



-m, 



/^y 



m.. 



m. 



\2m) (An) 2 \ fi 

Choosing m s = 550MeV and g = 13.5, the strength coefficients become C a = 96MeV and 
C s = 378MeV. 

In fig. 6 we show equipotential plots for the choices A = (graph a) and A = 1 (graph b), 
constructed by keeping two nucleons 1 fm apart and varying the position of the third one. 
As the plots are symmetric under rotation around the x-axes the specification of a single 
quadrant describes the spatial energy distribution. Inspecting this figure one learns that 
the predictions from both models are rather different, indicating that the effective seagull 
is very important. In the linear approach, the interaction is attractive over a wide region, 
whereas the non-linear scalar produces a repulsion around the triangular configuration and 
these differences should show up in observables. 

One is aware that a study based on just S trinucleon waves can provide only rough 
indications, since it is well known that D waves do play an important role in trinucleons. 
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Nevertheless, in the absence of a detailed study, we may assume that the trends associated 
with our equipotential plots would reflect in the binding energy. This assumption is sup- 
ported by the results of ref. |IB] where a three-body force given by our eq.(22) with A = 1 
was shown to produce a decrease of the binding energy which is rather welcome. 

As a final comment, it is important to point out that the discussion presented in section 
3 makes us to be biased towards the chiral scalar, but final conclusions must wait until a 
complete evaluation of the diagrams of fig 2c, which is now in progress. 
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Figure Captions 





Fig. 1. Interactions of nucleons (full lines), pions (broken lines) and scalar- isoscalar mesons 
(wavy lines) present in both linear and non-linear models (a) and only in the latter (b). 
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Fig. 2. Dynamical content of the effective scalar-isoscalar field in (a) nucleon-nucleon 
interactions, (b) the irNN kernel and (c) the seagull term; conventions are the same as 
in fig.l and the crosses in nucleon propagators indicate that they do not contain positive 
frequency components. 



18 



— <2> 



N ■ / \ / \ ✓ \ / 




/ \ II II II 

' * "I II II 

' V ' • i . . 



Fig. 3. The effective seagull (a) as composed by the irN — > %N (b) and irN — > ^iV (c) 
amplitudes. 
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Fig.4. (a) Amplitudes for the processes irN -> aN (A = 0) and tcN — ► SN (A = 1); (b) 
the corresponding irNN kernel. 
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Fig. 5. Basic diagram for the three-body force. 
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Fig. 6 Equipotential plots for the expectation value of the three body force in the trin- 
ucleon ground state, calculated using eq (36) with A = l.hGeV for the linear (a) and non 
linear (b) models. One of the nucleons is fixed at x — 0.5 fm, another at x — —0.5 fm and 
the position of the third one is varied. 
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